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ABSTRACT 

We present high sensitivity ^^CO and ^'^CO J=l— >0 molecular line maps covering the full extent 
of the parsec scale Haro 6-10 Herbig-Haro (HH) flow. We report the discovery of a molecular CO 
outflow along the axis of parsec-scale HH flow. Previous molecular studies missed the identification 
of the outflow probably due to their smaller mapping area and the confusing spectral features present 
towards the object. Our detailed molecular line study of the full 1.6 pc extent of the optical flow 
shows evidence for both blueshifted and redshifted gas set in motion by Haro 6-10 activity. The 
molecular outflow is centered at Haro 6-10, with redshifted gas being clumpy and directed towards 
the northeast, while blueshifted gas is in the southwest direction. The molecular gas terminates well 
within the cloud, short of the most distant HH objects of the optical flow. Contamination from an 
unrelated cloud along the same line of sight prevents a thorough study of the blueshifted outflow lobe 
and the mass distribution at the lowest velocities in both lobes. The cloud core in which Haro 6-10 is 
embedded is filamentary and flattened in the east-west direction. The total cloud mass is calculated 
from ^'^CO J=l-^0 to be ~ 200 Mq. The lower limit of the mass associated with the outflow is ~ 
0.25 Mq. 

Subject headings: ISM: clouds - ISM: jets and outflows - ISM:Herbig-Haro objects - individual Haro 
6-10 - stars: formation 



1. INTRODUCTION 

In the last ten years, wide-field CCD camera observa- 
tions have shown that it is quite common for outflows 
traced by Herbig-Haro (HH) objects to attain parsec 
scale dimensions, at leas t an order of m agnitude larger 
than previously thought (jReipurth et al.lflQQ'iT ). It is now 
becoming clear that optical HH objects (which are lit up 
from the cooling behind shock fronts within the flows), 
near- infrared H2 jets (tracing shocked regions), and the 
swept-up CO bipolar molecular o utflows are intimatel y 
linked to one another (see for e.g. iBachiller et al1ll995D . 
With the development of large-format heterodyne arrays 
at millimeter wavelengths and the availability of " on-the- 
fly" (OTF) mapping capabilities, millimeter mapping of 
the full extent of some of the known parsec-scale HH 
flows is possible with increased data sensitivity. This 
allows a more thorough study of the observational prop- 
erties of the outflows and provides better constraints on 
the theoretical models of molecular outflow entrainment 
mechanism. 

iDevine et al.l (|1999D discovered a giant 1.6 pc long HH 
flow, centered on Haro 6-10 at a position angle of 222° 
delineated by HH 410 and HH 411 at its edges, and HH 
412 and HH 184E along the HH flow axis. Haro 6-10 
is located in the L 1524 cloud within the B18 region of 
the Taurus Molecular Cloud Complex. It is a binary sys- 
tem composed of an optically visible southern component 
(Haro 6-lOS), and its infrared companio n (Haro 6-lON) . 
High-resolution 3.6 cm VLA map ( Reipurth et al.ll2004h 
strongly suggests that the optical component, Haro 6- 
lOS, is driving the current outflow activity. The same 
set of data suggests that there might be an additional 
stellar companion to Haro 6-lOS, more closely bound to 
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it than Haro 6- ION. 

There has been other evidence for the outflow activ- 
ity from the Haro 6-10 system such as 2.12 /im molec- 
ular hydrogen emission, ob served at the position o f 
the northern companion JCarr 1990; Herbst et al.lfl995h . 
iMovsessian fc Maga kian (1999) found a small HH jet as- 
sociated with Haro 6-10 itself at a position angle o f 195° 
with blueshifted radial velocities. iReipurth et al.l ()2004D 
reported a VLA jet along the same angle. These latter jet 
angles are significantly different from the 222° angle of 
the parsec scale HH flow. The difference in these position 
angles have led some authors to suggest that the outflow 
from Haro 6-10 has undergone precession o r reorientation 
(IReip urth et~al] 120041 : lOevine et al.lll999l ). iDevine et all 
(1999) have argued that HH 184F-G and HH 184A-B are 
two different flows originating at the same time from two 
different protostars, probably due to an interaction in 
the binary pair. They have position angles of 162° and 
231-249° respectively. HH 184E is the only knot in the 
HH 184 that is found along the parsec scale 222° axis 
and is the furthest away from driving source. 

In the past, a few attempts have been made to de- 
tect molecular CO outflow from the Haro 6-10 source. 
Small maps, only a couple of arcminutes wide, have 
been made toward Haro 6-10 which presented the ev- 
idence of weak non-Gaussian CO wings emiss ion in 
the d i rection of the o ptical object (Edwards fc~S"nel]| 
1981 iLevreaulf] [l988l : lHogerheiide"etaII fl997l . II998I : 



Chandler et al.lll998f) . Hitherto. Haro 6-10 has been con- 



sider ed to be devoid of molecular outflow (jPevine et al.l 
[1991). 

We made large, sensitive, ^^CO J=1^0 and ^^CO 
J^l^O maps of the Haro 6-10 region, over the full extent 
of the optical parsec-scale HH flow. Here we show the ev- 
idence for the redshifted high velocity CO gas along the 
uonatflaiii axis defined by the HH parsec scale flow. Due 
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Fig. 1.— The I^CO and ^^CO spectra averaged in 15' area centered on Haro 6-10. Foreground contamination from another cloud 
Taurus gives rise to the blueshifted emission peaking at ~ 4.5 kms~^. 



Fig. 2. — (See f2.jpg) Overview of the Haro 6-10 region. A zoomed view into the B18 region derived from a much larger-scale map from 
the Taurus molecular cloud survey ( NarayananiEin|2003) in ^^CO J=1^0 {a,b) and ^^CO J=1^0 {c,d). These images demonstrate the 
presence of two emission components at different velocities, along the same line of sight. Panels a and c show emission unrelated to the 
B 18 cloud, while panels b and d trace the B 18 cloud sub-structure. In panel 6 we identify three clouds found in B 18: L 1529, L 1524, 
B18w, from east to west). Several YSOs and T-Tauri stars are embedded in the L 1529 region. L 1524 hosts Haro 6-10 star and nearby 
HH 414 IRS (filled stars). B 18w at the western most point, hosts IRAS 04239-1-2436 source which drives parsec scale HH 300 flow. In the 
panel a we label Haro 6-10 binary system, the driving source of parsec-scale HH flow with terminating points at HH 410 toward SW and 
HH 411 toward NE. Contours at each panel start at 2a level and go in steps of Aa. 

to "contamination" by emission from another molecular 
cloud along the same line of sight, a thorough study of 
the blueshifted outflow lobe is limited. 



2. OBSERVATIONS 

A full mapping of ^^CO and ^■^CO in the J=1^0 tran- 
sition was performed with the SEQUOIA receiver at the 
Five College Radio Astronomy Observatory (FCRAO) 
14 m telescope during a period spanning 2003 and 2004. 
The receiver was configured as a dual-polarized 4x4 ar- 
ray. The orthogonal polarizations of the SEQUOIA ar- 
ray were averaged to produce spectra with higher signal- 
to-noise ratio (S/N). The telescope's half-power beam 
widths are 45" and 47" for ^^CO and ^^CO transitions 
respectively. Channel maps as well as individual spectra 
were checked for any scanning artifacts, baselined and re- 
gridded to the 22.5" sampled grid. RMS noise weighting 
was used to combine the data. For all calculations, the 
antenna temperatures were corrected for the main beam 
efficiencies of 0.45 for ^^CO and 0.5 for ^^CO. 

Both ^^CO and ^^CO spectra are smoothed to 
0.13 kms~^ channel spacing. The system tempera- 
tures (Tsys) in our observations range from 400 - 700 
K for 12 CO and between 200 - 500 K for ^^CO. Re- 



gions mapped with the higher noise level were repeated, 
combined and averaged in order to get a constant lower 
noise level over the whole extent of the map. The re- 
sulting mean rms per velocity channel is 0.2 K for ^^CO 
and 0.1 K for ^^CO. Antenna pointing and focus were 
checked every few hours and corrected using SiO masers. 

The analysis were done both using the Gildas software 
package and IDL software of Research Systems Inc. De- 
tailed studies of the physical parameters characterizing 
the outflows were performed using IDL. 

3. MORPHOLOGY OF THE CO GAS 

3.1. Two component emission in B18 cloud 

In Figure[T]we show ^^CO and ^^CO spectra averaged 
over a 15' area centered on Haro 6-10. In both lines, 
the averaged spectra consists of two Gaussian-like com- 
ponents: a stronger component centered at 6.4 kms^^ 
and a weaker component at '--^ 4 kms "^. A similar dou- 
ble pe ak profile has been observed bv lArce fc GoodmanI 
(|2001l hereafter AG) by studying the parsec scale out- 
flow from IRAS 04239-1-2426 in the adjacent B 18w cloud. 
They named the component at lower velocity "cloud A" , 
and surmised that is was produced by another cloud in 
Taurus along the same line of sight. 
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Fig. 3. — L 1524 cloud emission using our high sensitivity ^"^CO J=l— »0 data integrated in two velocity bins, representing two distinct 
component in L 1524 cloud. Contours at each panel start at 5cr level and go in steps of 10 cr in panel a (0.4 Kkms~^) and in steps of 20(t 
(0.9 Kkms~^) in panel b. HH 184 coordinate is taken from Bo Reipurth's General Catalog of HH Objects l |http: / / casa.colorado. edu /hhcat / |. 



The stronger Gaussian component is associated with 
the Haro 6-10 host cloud, L 1524. This is in agree- 
ment with the a mmonia core identif ied at LSR veloc- 
ity of 6.4 kms~^ (|Codella et al.lfT99l towards L 1524. 
Haro 6-10 is located at t he edge of the co mpact ammonia 
core (see Figure 3, iCodella et al.lll997l) . Our spatially 
averaged spectra towards the L 1524 cloud show some 
evidence for line wing emission at both blueshifted and 
redshifted velocities, although the contaminating cloud 
component at ^ 4 kms~^ obscures much of the wing 
emission at blueshifted velocities. Even with our higher 
sensitivity data, we do not see evidence for high velocity 
molecular gas as is seen in other prototypical outflows. 
Since the outflow is composed mainly of gas with slow ra- 
dial velocities, emission from outflowing gas at the slow- 
est velocities might be hidden under the ambient cloud 
emission. 

In order to examine this contaminating emission at 
4 kms^^, we used a 2° by 2° portion of the much larger 
FCRAO Taurus Molecular Cloud Survey in both ^^CO 
and "CO J=1^0 (Narava nan et al.ll2007f) . In Figured 
we plot for each line, the integrated intensity emission 
in two narrow velocity intervals. The first velocity in- 
terval preferentially selects the contaminating emission 
in the 3.5 to 4.5 kms""'^ velocity range, (Figure [Da, c) . 
This contaminating cloud component is traced by ^^CO 
molecule as an arc like filament, stretching from south of 
Haro 6-10 to north and appearing to bend toward IRS 
04239-1-2436 in B18w cloud and south from it again, (see 
Figure ^p) . The maximum intensity of this emission is 
south from IRS 04239-f 2436 source in the B 18w cloud, 
while at the location of Haro 6-10, the contaminating 
emission is quite clumpy. The ^^CO tracer shows sim- 
ilar morphology of the emitting gas in the region (see 
Figure [He) . 



In the other velocity interval, from 5.5 to 7.0 kms~^ 
we preferentially select the CO emission associated with 
the B 18 clouds, which have LSR velocities around 

6.4 kms-^ In ^^CO (see Figure Et), the B 18 cloud is 
revealed as a flattened, elongated structure in the east- 
west direction. In the same panel we identify L 1529, 
L 1524 and B 18w clouds within the B 18 region. 

The dramatic change in the morphology for each iso- 
tope at only a slightly different velocity offset indicates 
quite clearly that we are looking at two spatially distinct 
and possibly unrelated cloud features. From Figure [21 
it can be seen that towards Haro 6-10, the intensity of 
the contaminating emission is not as high as the inten- 
sity of the L 1524 cloud. However, towards B 18w, it 
is seen that the emission in both velocity bins seem to 
be of comparable intensity. The same was evident in 
the average spectra of the region presented by AG. In 
this study by AG, it was thus a clear cut in the case of 
IRAS 04239-1-2426 to concentrate only on the redshifted 
emission in the outflow. From the large-scale CO maps 
derived from the Taurus survey, it appears that while 
the blueshifted component in the Haro 6-10 spectra at 

4.5 kms~^ is part of the contaminating emission seen to- 
wards B 18w, the contamination is not nearly as severe 
as in the case of B 18w. Some of the emission towards 
Haro 6-10 at 4.5 kms""'^ might be part of the outflow 
system from that object. 

3.2. The Haro 6-10 cloud core 

In the following sections, we will examine in more de- 
tail the behavior of the ambient, blueshifted and red- 
shifted gas using our new ^^CO and ^^CO J=1^0 data 
of Haro 6-10 outflow, which are of smaller scale compared 
to the Taurus survey, but have better sensitivity. For the 
rest of this paper, the Taurus Survey data are not used 
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Fig. 4. — ^^CO line integrated in different velocity intervals. Upper row: Blueshifted gas. In panel a contours start and increase with 
2(7 certainty (2 Kkms~^), while in panel h they start and increase in 3(t steps (2.1 Kkms~^). Middle row: Line Core. Contours in c and 
d panels start at lOcr and go in step of 10 a (1 Kkms~^ at panel c and 1.2 Kkms~^ at panel d). In panel e contours start at 80" level 
(0.56 Kkms~^) and increase in the same step. Lower row: Redshifted gas. At both panels contours start and increase with 2a certainty, 
which is 0.2 Kkms^^ in panel / and 0.32 Kkms"'^ in panel g. In a, h we mark different knots of HH 184 IjPevine et al.in"999I V while in the 
rest of the panels we use HH 184 coordinate as defined in Figure |3] 
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Fig. 5. — (See f5.jpg) Best case scenario for the outflow in the Haro 6-10 region: blueshifted emission is integrated in the velocity interval 
from to 2 kms~^ and contours start at 0.3 Kkms~^ and go in steps of 0.1 Kkms~^; redshifted emission is integrated in velocity interval 
from 9 to 14 kms"'^ and contours start a t 0.5 Kkms~^ and go in steps of 0.5 Kkms~^. Molecular gas is overlayed on the optical image of 
the Haro 6-10 region IIDevine "etaLllTQQgl) A gain, we mark A, B, E, G knots of the HH 184 source IIDevine et al.|[T999l) . 



for any analysis. As already discussed, in the velocity 
range between 3.5 and 4.5 kms~^, the emission is domi- 
nated by the contaminating cloud component. Figure [3] 
shows the ^'^CO in two different velocity bins. Here, the 
contamination is clearly detected as strong feature at the 
northwest edge of the map (see Figure[3l3). In addition, a 
small clump of emission centered at Haro 6-10 is present, 
which is elongated and extended toward north, at the ~ 
30° angle with respect to the axis defined by the HH flow 
(crosses in Figures[3^). In the next panel, ^^CO emission 
traces the column density distribution within the L 1524 
cloud. The flattened structure of the L 1524 cloud seen 
in the larger-scale map of Figure [2] is also visible here, 
and Haro 6-10 lies close to the peak of integrated in- 
tensity emission. T his peak trace s the s ame core identi- 
fied in ammonia bv lCodella et al.l (|19971. The elongated 
L 1524 cloud is perpendicular to the HH outflow axis. 
From Haro 6-10, along the HH flow axis, the integrated 
intensity decreases, and HH 411 and HH 410 are found 
in the regions of very low ^'^CO integrated intensity. The 
ridge of ^'^CO emission is seen to extend northward from 
Haro 6-10 towards HH 414 IRS. 

3.3. Outflow Emission 

In Figure IH we show integrated intensity emission of 
^^CO in several different velocity bins. In the panels 
a, f and g we show the distribution of emission at the 
highest blueshifted and redshifted velocities. We believe 
that at these velocities we are selecting almost exclu- 
sively emission that arises from the outflow associated 
with Haro 6-10. Figure S] a shows that the blueshifted 
emission in the velocity range from to 2 kms^^ is spa- 
tially restricted and centered on Haro 6-10. Compared 
to the bubble-like structures defined by the knots A-G 
of HH 184, the blushifted emission has a smaller spatial 
extent and is elongated in the direction defined by A and 
F knots of HH 184. The spatial resolution of our data 
is such that is does not allow a more detailed compari- 
son. In Figure |4] /, 3 the redshifted emission consists of 
two bright regions centered on Haro 6-10 and HH 412 
respectively. The clump centered at Haro 6-10 has ex- 
tended emission toward southeast, similar to the highest 
velocity blueshifted gas identified in Figure [4] a. The 
redshifted emission seen toward HH 412, located about 
7' northeast of Haro 6-10, is detected to velocities as high 
as 14 kms~^, nearly 8 kms~^ from the ambient cloud line 
center. The emission is slightly extended along the HH 
flow axis. In the optical images (De vine et al.ill999i) . HH 
412 appears as a 1' long fairly diffuse emission feature, 
also elongated along the flow axis. Therefore, morpholo- 
gies of the CO and optical emissions are quite similar 
around HH 412. 

In Figure [5] w e show the opt ical Ha image of the 
Haro 6-10 region (jPevine et al.lfl99S) . overlayed with the 
^^CO J=l-^0 integrated emission in 9 to 14 kms~^ ve- 
locity range for the redshifted gas and to 2 kms~^ for 
the blueshifted gas. In the optical image, HH 410 and 
HH 411, lie at the extremes of the HH flow, and are seen 
against a background with many galaxies, demonstrat- 
ing that the flow has bursted out of the L 1524 cloud in 



both directions, and is possibly now moving through a 
region with greatly reduced column density. Our ^'^CO 
maps with negligible emission in these regions support 
this hypothesis. The expanded image of the CO outflow 
in Figure O shows that even toward Haro 6-10, the red- 
shifted and blueshifted outflowing gas are slightly offset, 
suggesting that the outffow is bipolar with the redshifted 
emission to th e northeast and bl ueshifted emission to 
the southwest. iDevine et al.l (|1999D speculated, based on 
morphology and brightness, that the large HH ffow has 
the same bipolar orientation. Additional information is 
provided by the radial velocity measurement of the small 
HH jet oriente d at position angle 195° that sh ows that it 
is blueshifted (jMovsessian fc Magakianlll999t ). and thus 
has the same velocity sense at the CO outflow, although 
the jet is oriented about 30° from the large HH flow. 

The distribution of CO emission in the near outflow 
wings, 2 to 3 kms^^ for the blueshifted gas and 8 to 

9 kms"^ for the redshifted gas, is shown in Figure[3] b,e. 
On the blueshifted side, emission is seen toward Haro 6- 

10 (as in Figure |4] a) along with much brighter and very 
extended emission oriented in a north-south direction. 
The southern most part of this emission lies on the HH 
flow axis and may be associated with the flow. However, 
most of the emission presumably has nothing to do with 
the flow. Therefore the emission in the 2 to 3 kms~^ 
range may be a mix of emission from the contaminating 
gas component as well as emission from the outflow. It is 
interesting to note that at these velocities the CO emis- 
sion close to Haro 6-10 is extended in the direction toward 
knots F and G in HH 184. Because of the contamination 
in this velocity interval we can not be sure whether this 
is outflow related or unrelated emission. Likewise in the 
8 to 9 kms~^ range there is emission near Haro 6-10 that 
may be related to the CO outflow. An additional feature 
is found to the north-east of Haro 6 -10 and conne c ts HH 
414 IRS to HH 413 and HH 412. iDevine et all (jHH) 
detected a small jet originating from HH 414 IRS (IRAS 
04264-1-2433) that was directed toward HH 413. They 
also suggested that the morphology of HH 413 knots re- 
sembled bow shocks originating from the direction of HH 
414 IRS. The morphology of the CO emission is sugges- 
tive that it is tracing the molecular counterpart to the 
HH 414 IRS flow. There could be also some contamina- 
tion at these velocities from the overall Haro 6-10 flow 
oriented north-east to south-west. It is noteworthy that 
the outflow axis, defined by the HH 414/413 optical flow 
and by the CO emission connecting them, is nearly per- 
pendicular to the axis of the main parsec-scale HH out- 
flow. Bright CO emission is also seen north-west of Haro 
6-10 that is presumably unrelated to either flow. Again 
the emission in this velocity range is likely a mix of out- 
flow emission and emission from the ambient gas. 

Finally it is interesting to note the structure in the 3 
to 5 kms~^ interval. Figure UJc. This velocity interval 
would be expected to include the outflow emission from 
Haro 6-10, since it lies a few kms^^ from the line center 
of the L1524 cloud emission, however the second cloud 
contaminates the emission. It is intriguing that there 
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Fig. 6.— l^co (thick) and I^CO (thin) spectra averaged in the 
redshifted lobe of the Haro 6-10 outflow. 

a finger- like emission from Haro 6-10 toward the north- 
east, and a bow shaped featm'e centered on Haro 6-10, 
and opening toward HH 410 in the southwest. These 
two features line up with the HH flow direction. The rel- 
evance of this structure to the outflow is unclear. In the 
following velocity bin (5 to 8 kms~^), the ^^CO emis- 
sion is dominated by L 1524 ambient cloud emission and 
no clear structure can be identified, although there is a 
maximum of emission towards Haro 6-10. 

4. MASS AND ENERGETICS 

The emission from ^^CO is usually optically thin and 
therefore a good tracer of the gas column density. How- 
ever, the emission can also be very weak, limiting the spa- 
tial and velocity extent to which outflows can be traced. 
While ^^CO is more readily detected, it is often opti- 
cally thick, even in the high velocity emission of out- 
flows. To derive accurate gas column densities we may 
need to correct for the optical depth of this line. Several 
recent studies have combined ^^CO and ^'^CO J=l— >0 
data to estimate the mass of t he outflow , using a velocity- 
dependent opacity correction (iBallv et al. 199 9l:lYu et al.l 
Il999t lArce fc GoodmanI 120011 : IStoiimirovic et al.ll2006D 
This approach usually requires construction of spatially 
averaged spectra, since ^'^CO emission is usually not de- 
tectable in most mapping positions. The application of 
this method in Haro 6-10 has some limitations due to the 
weak emission in the line wings. 

4.1. Outflow Mass 

In Figures[6]and[7]we show the spatially averaged ^^CO 
and ^'^CO spectra in the redshifted and blueshifted gas 
respectively. These spectra were constructed using the 
following procedure to isolate the strongest outflow emis- 
sion. For the blueshifted gas we found the integrated 
intensity at each location in our ^^CO map in the ve- 
locity range from to 2 kms^^. We then averaged the 
spectra at every location where the integrated intensity 
of "'^^CO was at least one-third of the peak integrated 
intensity. The ^'^CO spectrum was obtained by averag- 
ing the same spatial locations as was used to form the 
^^CO spectrum. For the redshifted gas we applied the 
same approach integrating over the velocity range of 9 
to 12 kms~^. The resulting spectra differ from the aver- 
age spectrum shown in Figure [T] in having both a weaker 
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Fig. 7.— i^co (thick) and ^^CO (thin) spectra averaged in the 
blueshifted lobe of the Haro 6-10 outflow. 

contaminating component and enhanced outflow wings. 
The outflow in ^^CO can be traced over a full velocity 
range from to 14 kms~^. 

Despite this averaging approach, the ^'^CO emission 
outside the line cores, has low signal to noise. Examining 
both Figures [6] and [7| we see that at the full spectral reso- 
lution, the ^'^CO emission is lost in noise at velocities be- 
low than 2 kms~^ and greater than 9 kms~^. The ^^CO 
outflow emission is too weak to use velocity-dependent 
opacity correction. We have instead integrated the emis- 
sion in four velocity intervals. In the near- wings of the 
redshifted outflow emission (between 8 and 9 kms~^) the 
ratio of ^^CO to ^^CO is 112 ± 20, consistent with the 
^^CO emission being optically thin. At higher outflow 
velocities (between 9 and 12 kms"'^), ^^CO emission is 
not detected, and the la lower limit on this ratio is 47. 
The ^'^CO emission in the blueshifted gas is even weaker, 
and only limits on the ratio can be obtained. In the near 
blueshifted outflow (2 to 3 kms~^) the 1 a lower limit is 
45 and at higher blueshifted velocities (0 to 2 kms~^) the 
lower limit is 20. Thus, where ^^CO emission is detected 
it is consistent with optically thin ^^CO emission and 
everywhere else we can only determine an upper limit 
to the ^^CO optical depth. Therefore to estimate the 
gas column density and outflow mass, we will assume 
that ^^CO emission is o ptically thin (equation 3 from 
nmirovic et al.l (|2006f) ). and' we will apply no correc- 
tion for the optical depth. 

The outflow mass as a function of velocity and position 
is then computed from M{v) = 2/imHANH2 (v), where 
fi = 1.36 is the mean atomic weight including He and 
other constituents, mn is the mass of the hydrogen atom 
and A is the physical area of one pixel at the distance of 
the source. Njjj is the molecular hydrogen column den- 
sity obtained u sing the relation N//^ = 1.1 x 10''Ni2 by 
iFrerking et al.l (1982) for the Taurus cloud; this result 
is consist ent with more r e cent d eterminations, summa- 
rized by iHarjunpaa et al.l ()2004f ) , for other nearby dark 
clouds. The greatest uncertainty in our mass estimates 
comes from the unc ertainty in the N/f, to N12 ratio 
(IFrerking et al.lll982D . 

The gas column density derived assuming that CO is in 
LTE depends on the gas temperature. For a gas temper- 
ature of 10 K we flnd 0.02 M0 in the velocity range from 
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TABLE 1 

Outflow Mass Estimates 





Tex = 10 K 


Tex = 25 K 




Mass 


Mass 


Velocity Interval 


Mq 


Mq 


Red 8 to 9 kms"! 


0.1 


0.15 


Red 9 to 14 kms'^ 


0.02 


0.03 


Blue to 2 kms-i 


0.01 


0.02 


Blue 2 to 3 kms-i 


0.02 


0.04 


Red Ibtal 


0.12 


0.18 


Blue Total 


0.03 


0.06 



9 to 14 kms^^, while the mass in the velocity range from 
to 2 kms^^ is found to be 0.01 Mq. For a gas temper- 
ature of 25 K, in the same velocity ranges, the redshifted 
gas has mass of 0.04 Mq and blueshifted of 0.02 Mq. 
We believe that these velocity ranges are dominated by 
outflow emission. 

The cutoff velocity of 9 kms~^ for the start of the 
redshifted component of the outflow and 2 kms~^ for 
the start of the blueshifted component correspond to ve- 
locity offsets from the systemic velocity of the cloud of 
2.5 and 4.5 kms~^, respectively. There is likely slower- 
moving gas at lower redshifted and blueshifted velocities 
that is missed in our analysis. Within the blueshifted 
velocities, the presence of the contaminating emission at 
4.5 kms~^ clearly makes it problematic to single out the 
low- velocity outflow component. Our outflow mass es- 
timates are therefore lower limits. In an effort to in- 
clude any lower- velocity outflow gas, we have repeated 
the mass calculation for the redshifted velocity range of 
8 to 9 kms~^ and derive a mass of 0.15 Mo, and for the 
blueshifted velocity range of 2 to 3 kms~^ and derive a 
mass of 0.04 Mq, both assuming a gas temperature of 25 
K. Our mass estimates are summarized in Table [1] 

The mass distribution with velocity for molecular out- 
flow has been shown t o have a power-law dependence, 
such that M(v) oc v"^ ()Richer et al.ll2000D . We have de- 
rived the mass- velocity relation for Haro 6-10 using the 
^^CO emission and assuming it is optically thin as we 
assumed earlier. We computed the mass in each velocity 
channel within the outflow and plotted this versus the 
velocity offset from the host cloud's mean velocity. The 
velocity range over which outflow gas is detected in Haro 
6-10 is very limited, making the determination of this re- 
lation very uncertain. In a log-log plot, the slope of the 
linear fit determines the 7 index. The value of 7 for the 
blueshifted gas is 6 ± 0.6 and for the redshifted gas is 
5.5 ± 0.6. 

4.2. Cloud Mass and Energy 

We determine the cloud mass by using the ^^CO map. 
Figure The line center optical depth at each point in 
the map is derived from the ^^CO peak temperature. The 
excitation temperature is obtained at each position by 
solving the radiative transfer equation for the excitation 
temperature, assuming ^^CO line to be optically thick. 
We find, the total mass of the cloud using this method 
to be ~ 200 Mq. The cloud mass estimate is mostly 
free from the foreground contamination, since we only 
searched for peak ^^CO emission in the narrow velocity 



range around line center. 

The kinetic energy of the cloud is dominated by the 
turbulent energy of the cloud, and is estimated using 
Etuib — 3/(16Zn2)McioudAv^. The mean turbulent ve- 
locity of the ambient gas is determined from the ^^CO 
line profiles in the cloud. We find the full line width at 
the half maximum to be Av = 1.2 kms^^. We find that 
kinetic energy of the L 1524 cloud is ^ 1.6 x 10^^ ergs. 

5. DISCUSSION 

Although at the lowest outflow velocities, the 
blueshifted emission is strongly contaminated by the 
foreground cloud emission, we see signatures of the out- 
flow from Haro 6-10 at the higher blueshifted velocities. 
There is significant overlap between blueshifted and red- 
shifted emission at the position of Haro 6-10, which to- 
gether with relatively low velocities observed in the line 
wings, suggests that the flow is in the plane of the sky. 
If the flow is in the plane of sky, slower radial velocities 
are hidden under the ambient cloud emission. In such a 
scenario, the clumps that we see defined in the redshifted 
gas (see for instance Figure [5]) are probably the highest 
velocity components seen projected in the map. 

A position- velocity (P-V) cut along the axis of the HH 
flow is shown in Figure [51 High velocity emission is seen 
at the location of Haro 6-10 (at offset in the P-V plot) 
and at the position of HH 412 (at offset 6.3). The highest 
velocity red-shifted emission is associated with HH 412. 
The P-V plot also reveals a general broadening of the 
CO emission at redshifted velocities between offsets of 
and -1-20', and at blueshifted velocities between off- 
sets of and —10'. The broadening at these positions is 
not symmetric about the cloud's systemic velocity, but 
rather, is wider at redshifted velocities for positive off- 
sets and at blueshifted velocities for negative offsets. The 
positive and negative positional offsets in Figure [8] corre- 
spond to the location of redshifted and blueshifted lobes 
respectively of the overall Haro 6-10 outflow system. The 
slight broadening at positive and negative offsets in the 
PV plot further bolsters the argument that the Haro 6- 
10 molecular outflow is, at least in the current epoch, 
a very low-velocity outflow that has slowed down close 
to ambient cloud velocities. If this is correct, it may be 
possible to see the subtle signature of the outflow at low 
velocities from constructing a centroid velocity map with 
velocities confined to the ambient cloud emission alone 
(6 to 7 kms^^). When such centroid velocity mapping 
was performed, for inst ance in the Cepheus A out flow 
system (see Figure 4 of iNaravanan fc Walkeil I1996D . it 
was found that the outflow signature is clearly seen as a 
velocity gradient along the outflow direction. We have 
performed such centroid velocity mapping analysis (fig- 
ure not shown) for Haro 6-10, but we do not see any such 
gradient present at ambient cloud velocities. This leads 
us to conclude that even if the outfiow has slowed to am- 
bient cloud velocities in Haro 6-10, the outflow has not 
imparted any significant velocity gradients over the nor- 
mally present random velocities inherent in the ambient 
cloud. 

Yet another puzzling feature of the Haro 6-10 outflow 
is, that it entrains very little gas. The morphology of the 
outflow is very clumpy and the total mass found in the 
outflow in only 0.25 Mq. In evolved sources, the previ- 
ous outflow episodes may clear up the ambient gas and 
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Velocity (km/s) 

Fig. 8. — Position- Velocity cut for the Haro 6-10 outflow, along 
the axis defined by the parsec-scale HH flow, position angle 40°. 
Haro 6-10 is located at an offset position of 0, and HH 412 is 
located at offset position of -1-6.3'. The cut is wide 5 pixels, which 
corresponds to '~2'. 

leave no gas to be entrained by later outflow episodes. In 
such a case, large ambient gas cavities are expected to be 
seen around source. Our ^'^CO maps show that no such 
cavities are found around Haro 6-10 and that the overall 
L 1524 cloud is somewhat flattened and filamentary, with 
the column density decreasing from Haro 6-10 source to- 
ward HH 410 and HH 411 along the outflow axis. While 
the L 1524 cloud is flattened, there appears to be enough 
ambient gas in the direction of the outflow lobes to be 
entrained (see Figure [S]). 

The clumpy morphology and low total mass may be a 
consequence of the outflow viewing geometry. If the out- 
flow is mostly in the plane of sky, much of its mass would 
be missed since we have excluded ambient gas compo- 
nents from our calculations. It would be interesting to 
obtain proper motion study of the HH knots to determine 
the orientation of the flow with respect to the plane of 
the sky. 

The outflow emission beyond Haro 6-10 is strictly con- 
fined to the HH 412 knot and it has the same elongated 
morphology as its optical counterpart. Therefore they 
seem to be created by the same mechanism. The HH 
knots are formed where the fast jet is interacting with 
the ambient gas. They can be found at the head of the 
jet or along its body at the locations of internal working 
surfaces where the new ejecta is catching up with the 
previously ejected slow-down gas. The CO clump at HH 
412 is detached from the Haro 6-10 source. It has small 
mass, and it is detected up to the projected velocities of 
8 kms~^ from the line center. 

Due to various contaminations, we had to use veloci- 
ties far from the line core to study the mass distribution 
in the case of Haro 6-10. This could be another reason 
why there is so little mass recovered in the outflow. Ob- 
servational studies show that outfiows show a power law 
distribution of the mass with respect to the velocity off- 
set from the cloud velocity. This distribution is such that 
most of the outfiow mass is at low velocities. The first 
mass point that we account for, in Haro 6-10 outflow, is 
offset by 2.6 kms^^ from the line center. For the Vlsb. 
of 6.4 kms~^ we calculated the outflow mass only for the 



points between 9-14 kms^^ and to 2 kms^^. Therefore 
we may be missing a significant fraction of the mass as it 
is increasing with the decreasing outfiow velocity, and the 
estimated outflow masses are consequently lower limit. 

Lastly, we can compare Haro 6-10 outflow to the HH 
300 parsec scale flow in the nearby B 18w cloud, studied 
by AG. AG found that the HH 300 outflow has a very 
clumpy structure and they identified five ^^CO redshifted 
clumps with masses from 0.03 to 0.23 Mq and radial ve- 
locities of about 3 kms^^ from the ambient cloud veloc- 
ity. They also noted that the clumps have different po- 
sition angles with respect to the driving source and con- 
cluded that the flow most likely precess with each clump 
corresponding to the different ejection event. The simi- 
larity of the HH 300 and Haro 6-10 outflow in terms of 
clumpy morphology and mass content is apparent. How- 
ever while the redshifted clumps of HH 300 flow are found 
at different position angles with respect to the IRAS 
source the redshifted clump at HH 412 location in Haro 
6-10 flow is found elongated along the jet (as is HH 412) 
axis defined by HH fiow. Although there are evidence 
that the current jet axis from the Haro 6 -10 differs from 
the a xis defined by the parsec scale fiow (iReipurth et al.l 
l2004HDevine et ai.lllQQalMovsessian fc Magakianlll999f) . 
the jet in Haro 6-10 must have had a stable orientation 
for a long time since all HH knots in the Haro 6-10 out- 
fiow which are found at large distance from the source 
(HH 410, 411, & 412) as weU as the HH 184E, lie on 
the 222° position angle jet axis. The HH 300 outflow 
exhibits larger impact on its host cloud, with ^■^CO line 
showing bipolar structure of the cloud gas at velocities 
close to the line center. Surprisingly, Haro 6-10 outflow 
does not seem to be affecting its host cloud's kinematics 
significantly. 

AG find a steep broken power law in the mass-velocity 
distribution with 7 — 4.0 at low outflow velocities and 7 
= 7.8 for outflow velocities greater than 1.85 kms~^. We 
find the slope 7 in the mass-velocity distribution of the 
Haro 6-10 outfiows to be around ~ 6. This value is much 
steeper than the averag e of 7 ~ 2 reported for a collection 
of outfiow sources (Ric her et al.l 2000). The steepening 
of the 7 index is expected for the older outflows, where 
once accelerated ambient material will slow down leading 
to accumulation of slow material. It is worth noting that 
the break in the mass- velocity distribution in the nearby 
HH 300 outflow happens very close to the line core. If 
the same holds for the Haro 6-10 we may be missing the 
shallower part at low velocities since we only account 
for the mass at velocities which have offset of at least 
2.5 kms~^ (toward redshifted) and 4.5 kms""'^ (toward 
blueshifted) from the systemic velocity of the cloud. The 
broken power law i s usually a good indicator of the jet 
entrainment model l| Zhang fc Zhengj[r997i ). 

6. SUMMARY AND CONCLUSIONS 

We made large, sensitive, ^^CO J=1^0 and ^^CO 
J=l— ^0 maps of the Haro 6-10 region, over the full extent 
of the optical parsec-scale HH flow. Here we summarize 
our main conclusions: 

1. The redshifted outflow component is clearly de- 
tected emerging from Haro 6-10 toward northeast, 
along the axis of optically defined parsec-scale HH 
flow. 
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Contamination from an unrelated foreground cloud 
along the same line of sight prevents a thorough 
study of the blueshifted outflow lobe, which we de- 
tect centered on Haro 6-10 in the opposite direction 
from redshifted lobe, and in the velocity range from 
to 2 kms-^ 

The mass and energies of the outflow are signif- 
icantly smaller then those of the host molecular 
cloud. However, most of the mass might be missed 
due to an unrelated contamination at lower outflow 
velocities, along the same line of sight, and/or in 
the plane of the sky orientation of the outflow. 



4. In the optical maps next to the HH 412, Haro 6-10 
parsec scale flow is crossed by HH 414/413 flow. In 
CO data at the lowest redshifted velocities, there 
seem to be evidence of HH 414 IRS activity toward 
HH 413. 
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